Stability of condensate in superconductors 
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According to the BCS theory the superconducting condensate develops in a single quantum mode 
and no Cooper pairs out of the condensate are assumed. Here we discuss a mechanism by which 
the successful mode inhibits condensation in neighboring modes and suppresses a creation of non- 
condensed Cooper pairs. It is shown that condensed and noncondensed Cooper pairs are separated 
by an energy gap which is smaller than the superconducting gap but large enough to prevent nu- 
cleation in all other modes and to eliminate effects of noncondensed Cooper pairs on properties of 
superconductors. Our result thus justifies basic assumptions of the BCS theory and confirms that 
the BCS condensate is stable with respect to two-particle excitations. 
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I. INTRODUCTION 

Seven decades ago London put forward the idea to ex- 
plain superconductivity by a rigid quantum wave func- 
tion covering all superconducting electrons. His bold vi- 
sion is included in all recent theories would the role of 
the wave function be fulfilled by the Ginzburg-Landau 
(GL) complex order parameter or the Bardeen-Cooper- 
Schrieffer (BCS) gap function 1 "^. The superconducting 
condensate is not called rigid but it is assumed to be 
sufficiently stable to form an effective vacuum of a new 
state, the state of broken symmetry 3 -. This stability is 
always used but its origin is not yet clear. 

In superconductors the macroscopic wave function is 
a phenomenological theoretical tool with a complicated 
underlying microscopic picture. In contrast, the macro- 
scopic wave function of dilute superfiuid gases is iden- 
tical to the intuitively clear Schrodinger wave function 
of the lowest energy single-particle state macroscopically 
occupied due to the Bose-Einstein condensation^ On su- 
pcrfluids we can outline the problem addressed here for 
superconductors . 

Landau has shown that the supercurrent is halted if ex- 
ternal perturbations excite bosons from the condensate 
to neighboring low-laying states. 5 According to this Lan- 
dau criterion a supercurrent in the ideal (noninteracting) 
Bose gas is unstable. Indeed, in the ideal Bose gas the 
energy of a neighbor state differs by the kinetic energy 
which is quadratic in momentum of this excitation. For 
any slow perturbation one then finds states of slower ve- 
locity vulnerable to excitation by Cherenkov-like mech- 
anism. In real Bose systems the situation is different. 
The interaction between bosons causes a reconstruction 
of the energy spectrum from the quadratic to the linear 
form of acoustic type. Perturbations slower than the cor- 
responding sound velocity then cannot excite particles so 
that the supercurrent is stable and the condensate wave 
function reveals the London rigidity. 



Excitations of the superconductors are principally dif- 
ferent. According to the BCS theory well supported by 
experimental experience, see e.g. Ref. 0, the supercon- 
ductivity is controlled by fermionic quasiparticles result- 
ing from broken Cooper pairs. We don't discuss here 
this familiar mechanism. Our central question is: Why 
Cooper pairs of nonzero momenta are not excited from 
the BCS condensate? 

Here we show that the multiple scattering corrections 
to the T-matrbs 7 - lead to gaps in the single-particle and 
two-particle energy spectra. The single-particle gap is 
the familiar BCS gap which is known to guarantee sta- 
bility with respect to the excitation of fermionic quasipar- 
ticles. We focus on the two-particle gap and show that it 
guarantees stability with respect to nucleation of super- 
conducting condensate in two or more momentum states 
and with respect to excitation of noncondensed Cooper 
pairs. These two gaps thus imply the stable condensate. 

The paper is organized as follows. In Sec. |H] we in- 
troduce the T-matrix approach first on a general level, 
Sec. Ill Al and then its simplified form for the separable 
interaction of the BCS type, Sec. Ill Bl Then we evaluate 
the T-matrix near the critical temperature for the con- 
densate in Sec. Ill C[ and for the noncondensed Cooper 
pairs in Sec. Ill Dl In Sec. IIIII we discuss the stability 
of the condensate with respect to nucleation of a second 
condensate, Sec. IIII Al and with respect to excitation of 
noncondensed Cooper pairs, Sec. IIII Bl Comparing the 
present approximation with the Kadanoff-Martin theory 
in Sec. IIII CI we point out the role of multiple scatter- 
ing correction in formation of the gap in the two-particle 
energy spectrum. Section IIVI is a summary. 



II. T-MATRIX APPROACH 

The noncondensed Cooper pairs are not covered by 
the BCS or Eliashberg theory. These approaches treat 
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Cooper pairs within the mean field which is nonzero 
only for states with macroscopic bosonic occupation - 
the states with the condensate. In our study we employ 
the T-matrix in which the pairing is independent of the 
condensation. We use the Galitskii-Feynman approxima- 
tion with multiple scattering corrections^ in modification 
of Ref . H 

A. General equations 

Let us introduce the theory. The full propagator is 
given by the Dyson equation 

Gf = + GrYl^-Gf, (1) 

where the selfenergy 

s t = -Jj-^2 a Qt ( 2 ) 

Q 

is a sum over four-momentum Q = (lu, Q) of interacting 
pairs, with Matsubara's frequencies u and discrete wave 
vectors Q corresponding to the sample volume CI. In the 
case of condensation mode Q is a four-momentum of a 
Cooper pair. <jq\ we call a Q-part. 

To avoid double-counts, the internal lines of the de- 
part of the selfenergy should not include processes related 
to the Q-mode itself. To this end we introduce the Q- 
reduced propagator 

G@i = G° + G°Y,qiGq_i, (3) 

which is dressed by all but the Q-part of the selfenergy 

Q'+Q 

The Q-part of the f selfenergy is obtained by closing 
the loop of \, line of the T-matrix by the reduced propa- 
gator 

a Qt (k)=T n (k,Q-k;k,Q-k)G m (Q-k). (5) 

The T-matrix is constructed from the Q-reduced prop- 
agators in the J, line and full propagators in the f line 

Tfi(k,Q-k;p,Q-p) = D(k,Q-k;p,Q-p) 

-^E^' ^- fc ; k',Q-k')G f (k')G m (Q-k') 

k' 

xT n (k',Q-k';p,Q-p). (6) 

Except for the reduced propagator, this is the standard 
ladder approximation. Here D can be either a general 
phonon propagator with vortices included or an effective 
interaction potential of the BCS type. The set of equa- 
tions is complete. 

The reduced propagator eliminates nonphysical re- 
peated collisions in the spirit of multiple scattering ex- 
pansion. If repeated collisions are not eliminated, which 



is achieved using approximation Gqi ~ Gj,, one recovers 
the original Galitskii-Feynman approximation. Impor- 
tance of the multiple scattering corrections can be seen 
from properties of the original Galitskii-Feynman approx- 
imation in the superconducting state. The T-matrix be- 
comes singular which signals the onset of pairing. The 
single-particle propagator, however, does not have the 
gap in the energy spectrumjSii£ With the multiple scat- 
tering corrections the gap develops^ 

The half-selfconsistent theory^ of Kadanoff and Mar- 
tin (KM) is recovered if we approximate the Q-reduced 
propagator by the bare one, Gq^ ~ G°. The KM the- 
ory yields the correct BCS gap but as noticed by Chen 
et aiii and confirmed below, the KM approximation re- 
sults in the ideal Bose gas of Cooper pairs. According 
to the Landau criterion the KM theory does not explain 
stability of the condensate with respect to excitation of 
noncondensed Cooper pairs. The reduced selfconsistency 
of the multiple scattering approach is thus essential for 
the excitation spectrum. 

B. Separable interaction 

For a discussion in this paper we employ the simple 
BCS interaction 

D{k,Q-k;p,Q-p) = -V6(uD-\e(k)\) 
x6(u;n-\e(Q-k)\)8(un-\e(p)\)0("v-nQ-P)\)-(7) 

From Eq. ^ one can see that this separable potential 
implies the separable T-matrix 

T n (k,Q-k;p,Q-p) = -T Q 9(Lo B -\e(k)\) 
xff(wD-|e(Q-*)|)ff(wD-|e(p)|)e(wD-|e(Q-p)|) (8) 

and Eq. ([6]) simplifies to a scalar equation 

4 = 7 + ¥"E G t(fc)G ,(Q-fc). (9) 

The sum over k is restricted by cutoffs of the BCS model. 
The Q-part of the selfenergy simplifies to 

VQtik) = -T Q G m (Q-k)e(u -\e(k)\) . (10) 

Equations Q-© and (|9|l-([T0)l form a closed set. 

C. Condensation mode 

The condensation of Cooper pairs happens in a single 
mode. Below we prove this assumption. From now on 
we reserve the index Q for the condensation mode, while 
the other modes will be denoted by Q 1 . The Matsubara 
frequency in Q is zero, Q — (0, Q). 

The T-matrix of the Q-mode diverges reaching values 
proportional to the volume f2, see Ref. [7]. To make a link 
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with the standard notation of the BCS theory we express 
this singular element as 



Tq = 



AA, 



(11) 



and split the selfenergy into the singular contribution of 
the Q-mode and the regular reminder 



E t (fc) = -AG Qi {-k)A 



a Q't( k ) 
- -AG Qi {-k)A + ^ m {k). 



(12) 



According to Eqs. ([T]) and the Q-reduced propagator 
relates to the full propagator as 



Gqi = Gi + GiAGq^-AGqi. 



(13) 



If one neglects renormalizations keeping only the gap, 
Ef « —AGqiA, i.e., E 0t w or G<^ « G°, this equa- 
tion becomes identical to the Nambu-Gor'kov equation 
with the BCS gap A. 

The T-matrix of the Q-mode 



1 



1 

V 



k B T 




Y^G t (k)G m (Q-k) 



(14) 



determines the gap. In the thermodynamical limit O — > 
oo, the T-matrix of the condensation mode diverges, i.e., 
I/Tq 0. Equation fll4j) then simplifies to the BCS-like 
gap equation 



° = 7 + ^rE G t(fc)G^(Q- 

fe 



k). 



(15) 



Gor'kov have analyzed equation (| 15[) close to the crit- 
ical temperature T c , where the gap is small. Keeping 
terms to the quadratic order in A he has shown that it 
leads to 



Ql 

2m* 



a + /3\A\ 2 = 0, 



(16) 



where m* is a Cooperon mass, /3 — 3/(2Ep) and a — 
-6Tr 2 klT c (T c -T)/(7( [3] E F ) are the GL parameters with 
n being the electron density, Ep the Fermi energy, and 
C[3] = 1.202 the Riemann zeta function. We restrict our 
attention to the vicinity of the critical temperature. The 
limiting form (|16[) will be thus sufficient for our discus- 
sion. 



D. Non-condensed Cooper pairs 

We expect that none of terms for Q' 7^ Q diverges with 
volume. This expectation is confirmed below. The Q'- 
reduced propagator then approaches the full one in the 
thermodynamical limit f2 — > oo, 



Gc 



Gf 



(17) 



and the T-matrix of a Q' ^ Q-mode 
1 1 k B T 



£G t (*)G+(Q'-A) 



(18) 



thus satisfies equation distinct from Eq. (|T4|. For the 
condensation mode (|14l) the gap enters only one of prop- 
agators while for noncondensation modes ([T8|) both prop- 
agators depend on the gap. This difference results in a 
suppressed excitation of noncondensed Copper pairs. 

The inverse T-matrix of noncondensation mode Q' = 
(0, Q') results from expansion of Eq. (TT5|) to the quadratic 
order in A as 



C 



Q' 2 

2m* 



2/?|A| 



(19) 



where C = 87r 2 fc|T c 2 /(7C[ 3 ]n). The factor of two in front 
of /3 follows from the fact that for noncondensed pairs 
both propagators depend on the gap. 



III. STABILITY OF THE CONDENSATE 

Now we are ready to solve the central problem of this 
paper. First we show that once the condensate is formed 
in the Q-mode, a parallel condensation in another Q'- 
mode is excluded. Second we show that the critical ve- 
locity for breaking the condensed Cooper into two quasi- 
particles is lower than the critical velocity of excitation 
of Cooper pairs out of condensate. 



A. Excluded parallel condensation 

The inverse T-matrix (fT9"| of the non-condensation Q'- 
mode cannot reach zero turning the Q'-mode into a par- 
allel condensation mode. To show this we first use the 
GL equation (I16j) , to express the T-matrix of nonconden- 
sation mode as 



C 
Tq> 



Q' 

2m* 



91 

m* 



(20) 



Values of the pair momentum Q are limited by the crit- 
ical current, Q 2 < Q 2 . The current is proportional to the 
square of the gap times the momentum, j cx Q|A| 2 . Us- 
ing Eq. (HU) one finds j cx Q (\a\ - Q 2 /2m*). The crit- 
ical current is the maximum one, <9j/9Q|q c = 0, which 
is achieved for Q 2 = 2m*\a\/3, see Tinkhami 2 -. Accord- 
ingly, 



Tq> 2m* 



Qr 



Q_ , H 

2m* 3 



(21) 



Inequality (|2"T1) implies that the mode of Q' ^ Q can- 
not become singular once the condensation develops in 
the mode Q. Therefore, a parallel condensation in two 
competitive modes is excluded. Briefly, there is a single 
condensate, as it is tacitly assumed in the BCS theory. 
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B. Excitation of Cooper pairs from the condensate 

Now we discuss a possibility to excite a Cooper pair 
out of condensate by an object moving with velocity v 
in the static condensate. Going into the floating coor- 
dinate system, this criterion is used to check stability of 
the condensate flowing with velocity — v around a static 
obstacle. 

The right hand side of Eq. (fT9"|) represents an energy 
of a noncondcnsed Cooper pair of momentum Q' . In the 
frame floating with the condensate, Q = 0, a Cooper pair 
can be excited from the condensate into a noncondensed 
state with the minimal energy cost |a|. Let us estimate 
under which conditions Cooper pairs can be excited by 
an external perturbation. 

According to the Landau criterion^ the external per- 
turbation moving with velocity v can excite the Cooper 
pair of momentum Q' if the Cherenkov condition 

Q 2 

V Q' = TT— + M (22) 
2m* 11 y ' 

is satisfied. This equation is solved by real Q' for 

M > M 

This velocity is higher than the critical velocity of pair 
breaking v c = A/fcp, where kp is the Fermi momen- 
tum. Indeed, from Eq. (fl"6|) follows A = y/\a\/P = 
yj |a|fcp/(3m), where m = m*/2 is the electronic mass, 
therefore 

|v| > VSv c . (24) 

Briefly, it is easier to break a Cooper pair into two quasi- 
particles than to excite it from the condensate into a 
noncondensed Cooper pair. 

C. Role of the multiple scattering corrections 

As pointed out above deriving Eqs. (fT4l and (TT51) . 
within the multiple scattering approach the propagators 
inside the T-matrix depend on the evaluated scattering 
process. As the condensation mode becomes singular, 
the two-particle propagations in the condensed and non- 
condensed modes become particularly different. Let us 



show that this difference is essential using the Kadanoff- 
Martin approximation which results from the present ap- 
proximation using G@i « G°, therefore it uses the same 
two-particle propagation Gf(k)G a (Q — k) for all modes. 

Now we confirm that the KM approximation results 
in the ideal gas of Cooper pairs^i. Since one of prop- 
agators is bare for all modes, the KM counterpart of 
equation ([19]) reads ^ = + H + /3|A| 2 so that 

equation (g0|) modifies to ^ = ^ - When the 

condensate moves, Q ^ 0, it is energetically favorable to 
start condensation in standing mode Q' = which stops 
the supercurrent. A similar problem appears for non- 
condensed Cooper pairs. In the frame moving with the 
condensate one finds the free-particle-like energy of non- 
condensed Cooper pairs Q' 2 /2m*. Therefore, according 
to the Landau criterion in the KM approximation the 
condensate is not stable. 



IV. SUMMARY 

We have discussed stability of supercurrents with re- 
spect to condensation in competitive modes and exci- 
tations of noncondensed Cooper pairs. It was shown 
within the Galitskii-Feynman approximation that mul- 
tiple scattering corrections yield the familiar BCS gap 
in the single-particle energy spectrum and also a smaller 
gap in the two-particle energy spectrum separating the 
noncondensed Cooper pairs from the condensate. This 
two-particle gap prevents parallel condensation of Cooper 
pairs in two or more modes. Moreover, due to the two- 
particle gap the critical velocity to excite noncondensed 
Cooper pairs is higher than the critical velocity of the pair 
breaking, therefore the noncondensed Cooper pairs do 
not affect stability of supercurrents. The present result 
justifies basic assumptions of the BCS theory in which 
the condensate is expected in a single mode and Cooper 
pairs out of the condensate are ignored. 
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